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1 Introduction 

In these notes, we will be mainly focussing on the proof of the so-called ASM-DPP conjec- 
ture of Mills, Robbins and Rumsey [22] which relates refined enumerations of Alternating 
Sign Matrices (ASM) and Descending Plane Partitions (DPP). 

ASMs were introduced by Mills, Robbins and Rumsey [21] in their study of Dodgsons 
condensation algorithm for the evaluation of determinants. DPPs were introduced by 
Andrews [1] while attempting to prove a conjectured formula for the generating function 
of cyclically symmetric plane partitions. 

1.1 ASMs from lambda-determinant 

The definition of the so-called lambda-determinant of Mills, Robbins and Rumsey [21] is 
based on the famous Dodgson condensation algorithm [12J for computing determinants, 
itself based on the Desnanot-Jacobi equation, a particular Plucker relation, relating minors 
of any square k + 1 x k + 1 matrix M: 

\M\ x IMJJ+Jl = iMjgl x \M\\ - |M* +1 | x \M£ +1 \ (1.1) 

where \M^f^'"f r \ stands for the determinant of the matrix obtained from M by deleting 
rows ii, i r and columns ji, j r . The relation (11.11) may be used as a recursion relation 
on the size of the matrix, allowing for efficiently compute its determinant. 

More formally, we may recast the algorithm using the so-called T-system (also 
known as discrete Hirota) relation: 

for any i, j, fceZ with fixed parity of i + j + k. Now let A = (aij)ij'e{i,2,...,Ti} be a fixed 
n x n matrix. Together with the initial data: 

T^mfi = 1 (£, m E 1*; £ + m = n mod 2) 

Ti,j,\ — a j-»+n+i i+j+n+i (i, j G Z; i + j = n + 1 mod2; \i\ + \j\ < n — 1) , (1.3) 
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the solution of the T-system (11.21) satisfies: 



T 0An = det(A) (1.4) 

Given a fixed formal parameter A, the lambda-determinant of the matrix A, denoted 
by \A\\ is simply defined as the solution To 5 o,n = |^4|a of the deformed T-system 

j-l,k + T i+ ij tk Ti_ij jk (1.5) 

subject to the initial condition (11. 3p . 

The discovery of Mills, Robbins and Rumsey is that the lambda-determinant is a 
homogeneous Laurent polynomial of the matrix entries of degree n, and that moreover the 
monomials in the expression are coded by n x n matrices B with entries bij e {0, 1, — 1}, 
characterized by the fact that their row and column sums are 1 and that the partial row 
and column sums are non-negative, namely 

k k 

^2bi d >Q >0 (k = 1,2,..., 7i- l;j = l,2,...,n) 

i=\ i=l 
n n 



i=l i=l 



Such matrices B are called alternating sign matrices (ASMs). These include the permu- 
tation matrices (the ASMs with no —1 entry). Here are the 7 ASMs of size 3: 




There is an explicit formula for the lambda-determinant [ 

nxnASMB i,j 

where Inv(-B) and N(B) denote respectively the inversion number and the number of 
entries —1 in B, with 



Inv(.B) = 6 <A* 




N(B) = 

Note that for A = — 1, only the ASMs with N(B) = contribute, i.e. the permutation 
matrices, for which lnv(B) coincides with the usual inversion number of the corresponding 
permutation, and therefore (I1.6P reduces to the usual formula for the determinant. 

Mills, Robbins and Rumsey [22] noticed that apart from the quantities hxv(B) and 
N(B), another "observable" of interest is the position of the unique 1 in the top row of 
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any ASM. We denote by t(B) the number of entries to the left of the 1 in the top row 
of B. 

Associating a weight 

W{B) = z t(B) y lnv(B)-N(B) x N(B) (L7) 

to each ASM B, we may form the partition function 

nxnASM B 

In the case n = 3 listed above, the ASMs receive respective weights: 1, zy, zxy, y, zy 2 , z 2 y 2 , z 2 y 3 , 
leading to Z^ SM (x, y, z) = 1 + zy + zxy + y + zy 2 + z 2 y 2 + z 2 y 3 . 

1.2 DPPs 

Descending plane partitions are arrays of positive integers of the form: 

0>1,2 0-1,3 ' ' ' ' ' ' °l,Mi-2 a l,Mi 

0-2,2 0-2,3 " " " " " " °2,^2 

03,3 ■ ■ ■ 0-3,^3 

Or,r ' ' ' O r ^ r 

such that the sequece /ij is strictly decreasing < fi iy and that for Aj = Hi — i + 1, 
A = 00: 

Oij > Oij + \ ciij > Oi + \j X{ < di t i < Aj_i 
for all i, j. By convention, the empty partition is a DPP. Here are the 7 DPPs of order 3: 

0, 2, 3, 3 1, 3 2, 3 3, 3 ^. 

The integers a^- are called parts. A DPP A is said to be of order n if a it j < n for 
all i, j. A part Ojj is said to be special if a^j < j — i. We denote by S(A) and NS(A) 
respectively the total number of special parts and the total number of non-special parts of 
any DPP A. Another observable of interest among the DPPs A of order n is the number 
of parts in A equal to the order, which we denote by M(A). To each DPP A of given 
order, we associate a weight 

W{A) = x SW y NS(A) z M(A) 

and define the partition function for DPPs of order n to be: 

Z<& p (x,y,z)= £ W(A) (1.10) 

DPP A of order n 

The 7 DPPs of order 3 listed above have respective weights: l,y, zy, zxy, zy 2 , z 2 y 2 , z 2 y 3 , 
as M(A) is the number of occurrences of the part 3, and the only special part is the entry 
1 in the fourth DPP. This leads to the partition function Z^p P = l + y + zy + zxy + zy 2 + 

z 2 y 2 + z 2 y 3 . 



(1.9) 
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1.3 The ASM-DPP conjecture 



The ASM-DPP conjecture as stated by Mills, Robbins and Rumsey [22j amounts to the 
identity between the partition functions of ASMs and DPPs as defined in the previous 
sections. This is the following: 

Theorem 1.1. The partition functions for the refined enumeration of ASMs and DPPs 
coincide, namely 

Z A sm(x, V, z) = Z DPP (x, y, z) 

This was finally proved in all its generality in jl] , and then generalized so as to include 
yet another observable in [5]. In the present note, we explain the rationale behind these 
proofs which strongly rely on manipulations of finite truncations of infinite matrices. 

For simplicity of exposition we shall start with the identity between the doubly refined 
partition functions Z^ M (x,y, 1) and Z ( ^ PP (x J y J 1). Each will be expressed as the deter- 
minant of a finite truncation (of size n x n) of an infinite matrix, and the identity between 
determinants will be derived from general principles relating the two "infinite" matrices. 
One key ingredient is the use of the double generating series for the matrix entries (see 
Appendix IA1 for definitions and properties). 

1.4 Outline 

The use of infinite matrices is somewhat non-standard in this context, and we would like 
to stress the power and beauty of the method. The infinite matrices occurring here actu- 
ally involve some fundamental object that came up in the study of so-called Lorentzian 
triangulations [8J, giving rise to one of the simplest examples of quantum integrable sys- 
tem. More precisely, random configurations of this particular class of triangulations may 
be generated by iterated powers of a transfer matrix T of infinite size. The problem was 
solved exactly by diagonalization of T in [8]. A drastic simplification of the problem 
comes from the existence of an infinite parametric family of such transfer matrices, which 
all commute with each other. 

The notes are organized as follows. 

In Section 2, we recall a number of facts about the transfer matrix of 1 + 1-dimensional 
Lorentzian triangulations, including other applications to trees and lattice path enumer- 
ation. 

In Section 3, we compute Z^ M (x, y, 1) by use of the Izergin-Korepin (IK) [T51 [T7] 
determinant formulation of the partition function of the bijectively equivalent configura- 
tions of the 6 Vertex (6V) model with Domain- Wall Boundary Conditions (DWBC). The 
difficulty here is to extract a homogeneous limit out of the IK determinant, and to put 
it in the form of the determinant of a finite truncation to size n x n of an infinite matrix 
which is independent of n. 

In Section 4, we compute Zp PP (x,y } 1) by use of the lattice path formulation of the 
problem [19], and by use of the Lindstrom Gessel-Viennot (LGV) determinant formula 
for the partition function of non-intersecting families of lattice paths. This expresses 
Z^ pp (x, y, 1) as the determinant of the finite truncation to size n x n of an infinite matrix 
independent of n. 
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In Section 5, we show how the relation between the double generating functions of the 
two infinite matrices above implies the identity between the determinants of any finite 
truncation thereof. This is the key to the proof of the ASM-DPP conjecture. We then 
show how this has to be adapted to include more refinements. 

In the conclusive Section 6, the ASM-DPP correspondence is placed in the wider con- 
text of the myriad of combinatorial objects and physical systems connected to ASMs. We 
also compare the two very different forms of quantum integrability underlying the ASM- 
DPP correspondence, one coming from the Lorentzian triangulations, the other from the 
6V model. 

We collect the useful formulas and definitions for generating functions and truncated 
determinants of infinite matrices in Appendix |A] 

Acknowledgments. These notes are largely based on work with E. Guitter, C. Krist- 
jansen, R. Behrend and P. Zinn- Justin. I thank the Mathematical Sciences Research 
Instiitute, Berkeley, California for hosting me while these notes were completed. 

2 The main actor: the transfer matrix of Lorentzian gravity 
2.1 1+1D Lorentzian gravity 

Discrete models for 1+1D Lorentzian gravity are defined as follows. They are statistical 
models whose configurations are discrete space-times, in the form of random triangulations 
with a regular discrete time direction (an integer segment £2] C N) and a random space 
direction, modeled by random triangulations of the unit time strips [t, t + 1], t G N, by 
arbitrary but finite numbers of triangles with one edge along the time line t (resp. t + 1) 
and the opposite vertex on the time line t + 1 (resp. t). All other edges are then glued 
to their neighbors so as to form a triangulation. Each "horizontal edge" along a time line 
t is shared by two triangles, one in each time slice [t — l,t] and t,t+ 1]. The boundary 
conditions along the time lines may be taken free, periodic or staircase-like [8]. A typical 
free boundary Lorentzian triangulation G in 1+1D reads as follows: 




> 

space 

These triangulations are best described in the dual picture by considering triangles as 
vertical half-edges and pairs of triangles that share a time-like (horizontal) edge as vertical 
edges between two consecutive time-slices. We may now concentrate on the transition 
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between two consecutive time-slices which typically reads as follows: 



(2.1) 



with say i half-edges on the bottom and j on the top (here for instance we have i = 9 
and j = 10). Denoting by \i) and with_j i, j G Z+, the bottom and top Hilbert space 
state bases, we may describe the generation of a triangulation by the iterated action of a 
transfer operator T with matrix elements ly = ( l+ i :> ) between states \i) and \ j). Note that 
the corresponding matrix T = (Tjj)jj e z + is infinite. We shall deal with such matrices 
in the following. As detailed in Appendix |A} a compact characterization of the infinite 
matrix T is via its double generating function: 



f T (u,v) = T i,j 
i,j>0 



1 



1 — It 



2.2 Integrability 

To make the model more realistic, we may include both area and curvature-dependent 
terms, by introducing Botlzmann weights w(Q) equal to the product of local weights of 
the form g per triangle (area term) and a per pair of consecutive triangles in a time-slice 
pointing in the same direction, either both up or both down (curvature term). The rules 
in the dual picture follows: 



g g a a 

For instance, in the example (12. ip above with i — 9 and j = 10, the product of local 
weights is g 19 a 9 . For the staircase-like boundary conditions of [Sj, namely assuming that 
each state-to-state transition as in (12 .ip has at least one leftmost half-edge on the bottom 
and one rightmost half-edge on top (and not counting the leftmost and rightmost half-edge 
weights g 2 ), it is easy to compute the new transfer operator 7(g, a) with matrix elements 

T(g,a) hJ = (ag)^ £ (* J a" 2 * (2.2) 



k=0 



for i,j>0 which expresses the transition between states \i + 1) and \ j + 1). Equivalently, 
the double generating function reads: 

fT( g ,a) (u, v)=J2 T(g, a\ jU l v ] = — — - — — (2.3) 

f-^ Q 1 — ga{u + v j — g z (l — a z )uv 

and we have T(g, a) = T(ga, ga, g 2 (l — a 2 )) in the notations of Appendix [A] 

1 Here and throughout the paper, we use the notation = {0, 1, 2, ...}. 
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This model turns out to provide one of the simplest examples of quantum integrable 
system, with an infinite family of commuting transfer matrices. Indeed, we have: 

Theorem 2.1. [8] The transfer matrices T(g, a) and T(g', a') commute if and only if the 
parameters (g,a,g',a r ) are such that <p(g,a) = <p(g',a') where: 

tp(g,a) = 4- } - (2.4) 

ug 

This is easily proved by using the generating functions. This case corresponds to the 
family of matrices T Stt (a) of Appendix lAl with s — 1, t — <p(g, a) and a = ga. 

2.3 Diagonalization 

To diagonalize T, we may first consider finite size truncations T' 0,fc l = (Tij(g, a))^ 
It is easy to see that the (k + 1) x [k + 1) symmetric matrix is diagonalizable, with 
eigenvalues Xf' k \g,a) = g 2l (l + 0(g 2 )), % G [0,k], all with formal series expansions in 
powers of g with coefficients in Z[a]. As k increases, we get more and more eigenvalues, 
with power series expansions that stabilize. In this sense, the limiting infinite matrix has 
an infinite set of eigenvalues Aj = g 2t (l + 0(g 2 )), i G Z + , with well-defined formal power 
series expansions in g. 

Setting (p(g, a) = q + q^ 1 for some q G C*, and introducing a new variable 

1 - q- x ga 
1 - qga 

we may rewrite the double generating function of T(g, a) as: 

1-Ag 2 



fT(g,a)(u,v) 



[1 — qu){l — qv) — \{u — q){v — q) 



oo 



V / W F (<7 -u) m (1- \q 2 \ y/T^g - v) m ^ ^ 



where we identify 



^ Q (1 - qu) m+1 \l-q 2 ) (1 - qv) m+1 



1 — \a 2 
1-q 2 



as the m-th eigenvalue of T(g, a) , m G Z + and 



f _ \r ,,N,,i _ yi-g 2 (g-«) m 

Jv(m)( u )- u i a - (\-qu) m + l 

as the generating function for the corresponding eigenvector Note that (v ( - m ^) m e% + 

form an orthonormal basis of the Hilbert space of states w.r.t. the standard scalar product 
u ■ v = J2iez + u i v i- It is eas Y to show that A*™) = g 2m (l + 0(g 2 )) as a formal power series 
of g, thereby proving that these are the limits of the eigenvalues of the truncated matrices 
as the size k —¥ oo. 

This was extensively used in [8] to compute correlation functions of top/bottom bound- 
ary loops in random Lorentzian triangulations. We want to stress here the very simple 
form of the generating function (12. 3ft . which will reappear later in these notes. 
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2.4 Trees 



For staircase boundary conditions, the dual random Lorentzian triangulations introduced 
above may be viewed as random plane trees. This is easily realized by gluing all the 
bottom vertices of parallel vertical edges whose both top and bottom halves contribute 
to the curvature term (no interlacing with the neighboring time slices). A typical such 
example reads: 




Note that the tree is naturally rooted at its bottom vertex. 

To summarize, we have unearthed some integrable structure attached naturally to 
certain plane trees. Note that in tree language the weights are respectively g 2 per edge 
(except the left/right boundary ones), and a per pair of consecutive descendent half-edges 
(from left to right) and per leaf. 



2.5 Paths 

There is yet another interpretation of the transfer matrix T(g, a) of 1 + ID Lorentzian 
gravity, in terms of lattice paths. First notice that T(g,a) = V (g , a)V t (g , a) for some 
(infinite) lower triangular matrix V(g, a) with entries 

V{g,a) i)k = g i a i - k r') {i,keZ + ) (2.6) 



and double generating function: 



fv(g,a)(u,v) = (2.7) 

1 — agu — guv 

In the notations of Appendix IA.21 we have V(g,a) = L(a _1 ,5 | a). 

Consider paths on the positive quadrant of the two-dimensional square lattice Z+ 2 , 
with steps (—1,0) (horizontal) and (0,1) (vertical), as illustrated in FigfjTa). Then the 
total number of paths from the point (i, 0) to the point (k, k) on the diagonal is Q) . 
Moreover if we attach a weight ga per horizontal step and g per vertical one, we get a 
total contribution of {ga) k ~ l g k (lj = V(g, a)^, which we interpret as the partition function 
for weighted paths from (i, 0) to (k, k). This quantity will also reappear later. 

Note that in this language T(g, a)ij is the partition function of lattice paths from (i, 0) 
to (0, j) in the positive quadrant, with weights ga (resp. g) per horizontal step below 
(resp. above) the diagonal and g (resp. ga) per step above (resp. below) the diagonal 
(see Figfllh) for an illustration). 



8 



(i,0) 



(i,0) 



(a) 



(b) 



Figure 1: A typical path contributing to the matrix element V(g,a)i i k (a), and to T(g,a)ij (b). These paths are 
taken on Z+, with steps ( — 1,0) and (0, 1). In the latter case, we have indicated the weights of the steps below 
and above the diagonal. 

This formulation allows to visualize immediately the truncated transfer matrix T^- ' k ^(g, a) 
as corresponding to paths within the square [0,k] x [0,k] C Z+. For such paths, both 
the portion below the diagonal and that above are within the same square, so that we 
may write T^°' k ^(g, a) = V^°' k ^(g, a)V^°' k ^(g, a)*, wich immediately yields the determinant 
of T^°> k \ as V^ ' k ^(g,a) is lower triangular (see also Appendix IA.4j) : 



This is compatible with eigenvalues A| ' (g, a) = g 2l (l + 0(g 2 )) for i — 0, 1, k. 

3 Enumerating ASMs 

3.1 ASMs, 6V model and the IK deteminant 

As discovered by Kuperberg [18] . ASMs of size n x n are in bijection with the so-called 
Domain- Wall Boundary Condition Six Vertex (6V-DWBC) model on a square grid of size 
n x n. The latter configurations are choices of orientations of the edges of a n x n grid 
of the (two-dimensional) square lattice, in such a way that at each vertex exactly two 
edges point to- and two point from- the vertex. Moreover oriented external horizontal 
(resp. vertical) edges are attached to the boundary vertices, in such a way that external 
horizontal edges point towards the grid and vertical ones from the grid. We display below 
the 6 possible vertex configurations ai, a2, £>i, 62, c±, C2 obeying the above rules, as well as 



det {T^%,a))=det (V^(g,a)) 2 = g 



,k(k+l) 
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a sample grid showing the external edge boundary condition: 



6V 



b. 



ASM ooo 



1 -1 



6V 
DWBC 




We have also indicated the dictionary between the vertex configurations and the ASM 
entries. It is easy to understand the bijection as follows. The ±1 entries correspond 
to vertices where both the horizontal and vertical flows (indicated by the direction of 
the edges) are reflected. The alternation of 1 and —1 entries corresponds to odd and 
even order flips along each row or column of the grid. The DWBC boundary condition 
ensures that the first and last encountered non-zero elements in the ASM along rows and 
columns must be 1, as one needs an odd total number of flips to reverse the external edge 
orientation, when going along a row or a column. 

The 6V model had been extensively studied in the physics literature. With a suitable 
parameterization of the Boltzmann weights Oj, 6j, Cj, the model forms the archetypical 
example of an integrable lattice model, as it admits an infinite family of commuting 
transfer matrices, that can be diagonalized for various types of bopundary conditions 
using the Bethe Ansatz techniques. These integrable weights are defined as follows. Each 
row (resp. column) of the grid carries a complex number Z{, i = 1,2, ...n (resp. Wi, 
i = 1,2, ...n ) called spectral parameter. Moreover the weights depend on a "quantum" 
parameter q E C*. We have the following parametrization of the weights: 

b(z, w) 



a{z, w ) 



qz — q w 



q z — qw 



c\z,w) 



(q 2 - t 2 : 



where a(z,w) is the weight for a vertex of type a± or a2 at the intersection of a line with 
parameter z and column with parameter w, etc. With this parameterization, the model 
has an infinite family of commuting row-to-row transfer matrices, and can be exactly 
solved by Bethe Ansatz techniques. Using recursion relations of Korepin pTj , Izergin [15] 
obtained a compact determinantal formula for the partition function of this 6V-DWBC 
model, defined as the sum over edge configurations of the product of local vertex weights, 
divided by the normalization factor YYi=i c i z ii w i) (to make the answer polynomial in the 
z's and w's). It reads: 



H^aiz^Wj^z^Wj 



det 



1 



where A(z) = H 1<{ 



<j<?i\ 



A(z)A(w) i<i,i<n \a(z i} Wj)b(zi } Wj) 
stands for the Vandermonde determinant of the z's. 



(3.1) 



3.2 Homogeneous limit and computation of Z^g M (x,y,l) 

In the above bijection between 6V configurations and ASMs, it is easy to track both 
quantities N(B) and lnv(B) in terms of 6V weights. We find that 

N — n N 
N(B) = N C2 = lmr(B) - N(B) = N ai = N a2 = 
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where N av N^, N c . stand for the total numbers of vertex configurations of each type. The 
determinant result above can therefore be used to compute the refined partition function 
Z^g M (x, y, 1) for ASMs, which counts ASMs with a weight x/y per entry —1 and a weight 
y for each inversion. Setting 

HD 2 »=(?)' <-) 

we have: 

Zf SM {x lVl l)= x N ^W nv{B) - N{b) =b- n ^Z^(a,b,c) (3.3) 

ASM B 

where Zgy (a, b, c) refers to the homogeneous limit of the partition function (13. ip of the 6V 
model in which all a(zi,Wj) tend to a, etc. This is obtained by letting all Zi — > r and all 
Wi — > r _1 , with a = a(r, r _1 ), b = b(r, r _1 ) and c = c(r, r" 1 ). This and more refinements 
were worked out in [4]. We have the following remarkable result: 

Theorem 3.1. The partition function for refined ASMs reads: 

zf SM {x,y,l)= det ({l-u)I + vG) (3.4) 

0<« ,j<n— 1 

where v is any solution to the equation 

xv(l-v)=v + y(l-v) (3.5) 

and the n x n determinant is the principal minor for the n first rows and columns of the 
infinite matrix Masm = (1 ~ ^)I + vG whose entries are generated by 

fM ASM ( u , v ) = T — — + i — t ^ — ( 3 - 6 ) 

1 — uv 1 — xu — v — [y — x)uv 

Note the remarkable similarity between the generating function for the matrix elements 
of G and the that of the transfer matrix for 1 + ID Lorentzian triangulations f 12 . 3 j) . The 
two actually match up to a rescaling u — > uj\fx and v — > v^fx (which amounts to a 
conjugation by the diagonal matrix y/xT) and upon identifying x = g 2 a 2 and y = g 2 . 
Note also that G = T(x, l,y — x) in the notations of Appendix IA. 21 

Let us now give a sketch of the proof of Theorem 13 . 11 The determinant (13. ip is 
singular in the homogeneous limit, but we may Taylor-expand the matrix entries within 
the determinant around the homogeneous point. For a = a(r, r _1 ), b = b(r,r~ l ) and 
c = c(r, r _1 ) this reads: 



r(n) u _ (ab) n J ^ (\ (I d* \ (I & \ c{u-\v) 



Z 6V '(a,b,c) = — — det 



c n o<i,j<n-i \ [ \i\ du l J \j\ dvi J a(u 1 ,v)b(u 1 ,v) 

Noting further that 



c(u 1 ,v) 



a(u l ,v)b{u 1 ,v) uv — q 2 uv — q 2 
and introducing the infinite matrices A± with elements: 

(! - jeZ+) (3 - 7) 

we have the following straightforward: 
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Lemma 3.2. We have 



A ± = r -2 ]_ q ±2 (U(a±,P±Y) 1 U(a' ± ,P' ± ) (3i 



for U(a,f3) the infinite upper triangular matrix with entries U(a,(3)i t j = (ya 1 ^ (see 
Appendix \A.ty) , and where the parameters read: 



1 — q 2 r 2 a q 2 — q 2 

~2 _ n 2 



1 ' n H H I 2 2 n nl 

«+ = , P+ = — T > a + = ~Q T P+ > & 



r* — q* a+ 

and the parameters with — index are obtained form those with + by the substitution 
q —> q^ 1 . 

Proof. The statement of the lemma is an immediate consequence of the fact that the 
Taylor-expansion expression (13.71) around (u,v) = (r -1 ,?" -1 ) turns into the following dou- 
ble generating functions for the matrix elements of A±: 



f A± (u,v) = (Atk 



1 



(r 1 + u)(r 1 + v ) — q ±2 
Moreover, using the generating function for the matrix elements of U(a, (3) of Appendix 

/[,W)(v)= i-Mi + H' 

and for a(3 ^ 0, U(a, = U(—l/{3, —1/a), we finally compute by convolution product: 

A/ t (a,/3)- 1 ^(a',/3')(' U ' t ') = fv(-l/$ ,-l/a)(v,u) * fu(a',/3')(u,v) 

dt 1 1 



2int 1 + - it- 1 ) 1 - P'v(l + aft) 
1 



The lemma follows from comparing this with the expressions for f A± (u,v). □ 

The Lemma is easily extended to finite truncations of the infinite matrices A± and 
U, as U is upper triangular and (t/*) -1 is lower triangular. Therefore we may write 
A± = r - 2 ^_ q ±2 (U± t ~-j-ijj\^' n !] -with the obvious shorthand notations. Going back 
to our original determinant, we find that 

Z${a,b,c) = @£det (aL '- 11 - 4°' n " 1] ) 
= M!! d et(A [Q ' n - 1] )det fl- r " 2 ~ g'V '"- 11 ^C/i '"- 11 yi^n-HV^-Hy 1 

where I stands for the n x n identity matrix. The last product of 4 matrices is finally 
identified with the matrix G of Theorem 13.11 and computed by use of generating functions, 
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non-special 




non-special 

8 8 7 7 6 
7 6 6 4 
4 4 
2 



special 

4 3 
1 



1 2 3 4 5 6 

Figure 2: Non-intersecting lattice path configuration for a sample DPP of order n > 8. We have indicated the 
domains in which horizontal steps correspond to special/non-special parts. 



whereas the determinant of ^l '™ -1 ' follows from its expression as a product of triangular 
matrices. Collecting all the factors finally yields fl3.4H3.6p . with the following identification 
of parameters: 



x 



Q Q 



-2 



q 1 r 



qr 



qr 



q 1 r 1 



q 1 r — qr 1 



v 



q A — q 



v 



g z — q~ 



4 Enumerating DPPs 

4.1 Lattice path formulation of DPPs 

The DPPs are in bijection with configurations of non-intersecting lattice paths illustrated 
in Figj2] and defined as follows. Like in Sect. I2.5[ the paths take place in the positive 
quadrant (Z + ) 2 , with the same steps but different weights and boundary conditions. The 
paths start along the x axis at positions of the form (sj, 0) (i = 1,2, ...,r recorded from 
right to left) and end along the y axis at positions (0, + 2) (z = 1, 2, r recorded from 
top to bottom). We add a final horizontal left step at the end of each path. Reading 
paths from left to right and top to bottom, we record the vertical positions y = a^j of the 
j-th horizontal step from the left taken on the z-th path from top (steps with y = are 
not recorded). These form a DPP with r rows, of order any n > s\ + 2. Conversely to 
each DPP with r rows we may associate such a path configuration. Note that the starting 
points are such that Sj = Aj — 1, where Aj = yUj — % + 1 the total number of parts in the 
row i. 

The special parts correspond to horizontal steps taken in the strict upper octant y > 
x + 1 of the plane, and the remaining parts correspond to the horizontal steps in the 
domain 1 < y < x + 1, while horizontal steps along the x axis do not count (weight 1). 
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4.2 Computation of Zp PP (x, y, 1) 

The computation of Z^ PP (x, y, 1) uses the Lindstrom-Gessel-Viennot [211 [H] determinant 
formula expressing the partition function for non-intersecting lattice paths with fixed 
atarting points and endpoints Z as a determinant det(Zjj) where Zjj is the partition 
function for a single path from the i-th starting point to the j-th endpoint. This leads to 
the following: 



Theorem 4.1. ^ The partition function Zp PP (x, y) for DPPs of order n with weight 
per special part and y per other part reads: 

Z ( D n) PP (x,y) = det{l + H^) (4.1) 

where the determinant is that of the finite truncation to the n first rows and columns of 
the infinite matrix (Mopp)i,j = Si,j + Hij, i,j G Z +; with generating function: 



x 



f MDP M v) = Yl (MnppkjuV = — L_ + - — (4.2) 

' 1 — uv 1 — u 1 — xu — v — \y — x)uv 



Again, note the close similarity between the matrix H and the transfer matrix T 
for 1 + ID Lorentzian triangulations. Our proof of the ASM-DPP conjecture will be 
based on this similarity. Note also that in the notations of Appendix IA.21 we have 
H = yS(J.-S)- 1 T(x,l,y-x). 

Let us now sketch the proof of Theorem 14.11 The sought after partition function is 
a sum over all configurations of r non- intersecting paths (0 < r < n — 1) with fixed 
r starting points (sj,0), i = l,...,r and endpoints (0, Sj + 2), i = l,...,r. According 
to the Lindstrom-Gessel-Viennot theorem, this is the sum over minors of the 

n x n matrix D whose entries Z\j is the partition function for a single path from (i, 0) to 
(0, j + 2). It also has the simple expression: 

n-l 

E E \D\Z::Z = det(i+D) 

r=0 0<s 1 <...<s r <n-l 

Such a path is split into three pieces: (i) between the x = axis and the first hit on 
the x = 1 axis (ii) between the x = 1 axis and the diagonal line y = x + 1 (iii) between 
the diagonal line y — x + 1 and the vertical axis y = 0. Each piece receives a specific 
weight, with total contribution: 

» Min(k,j+l) / . /-li- 
ft \ k _i I J + i 



fc=0 fcO 

where we have first summed over (i) paths from (i, 0)to (k, 1) with i — k horizontal steps 
along the x = axis and one final vertical step (ii) paths from (k, 1) to (£,£ + 1) on the 
diagonal y — x + 1, for which k — I horizontal steps must be chosen among a total of 
k, each weighted by x as these correspond to special parts (iiii) paths from (£,£+ 1) to 
(0, j + 2) for which £ horizontal steps must be chosen among a total j + 1, each weighted 
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by y as they correspond to non-special parts with one extra y factor for the additional 
final horizontal step. 

Note that the matrix elements of D are independent of n. We may therefore consider 
the extension of D to an infinite matrix D with matrix elements given by Djj of eqn. ()4.3p . 
for i,j G Z + . The theorem follows by identifying the infinite matrix H with D, and 
therefore H^ n ~^ with D. 

5 Proof of the ASM-DPP conjecture 

5.1 Proof of Z { ^ M (x, y, 1) = Z% PP (x, y, 1) 

The expressions H3.4[) and (14. ip for respectively the partition functions Z^g M (x, y, 1) and 
Z^ pp (x J y,l) are determinants of the principal minor of size n of some infinite matrix, 
in other words, these are the determinants of a finite truncation to the n first rows and 
columns of infinite matrices. 

There is a very simple relation (independent of n) between the generating functions of 
the two infinite matrices Masm and Mppp, namely: 

u 

(1 - ^(i-fl/MAsaK^ = (1 -U)(l- (1- v)v)f MDpp {u,v) (5.1) 

direct consequence of (13. 5ft . 
Let us translate this back into a finite matrix relation upon truncation. First, for 
any matrix A with generating function /^(m,^) the function fM(u,v) = (1 — au){\ — 
bv)fA(u,v) is actually the generating function of the infinite matrix M — (I — aS)A(I — 
bS t ) where S is the strictly lower triangular shift matrix with elements Sij = for 
i,j G Z + , and S l its strictly upper triangular transpose. Upon truncation to indices in 
[0, n — 1], we have the obvious relation (see Lemma [A. 2 1 in Appendix |A|) : M^°' n ^ 1 ^ = (I — 
aSy°' n ~ 1 ^ A^-°' n ^ 1 ^ (I — bS t y°' n ~ 1 \ due to lower triangularity of I— aS and upper triangularity 
of I — bS l . Note that both matrix truncations are unitriangular, hence have determinant 
1 so that det(M[ '"- 1 ]) = det^ '™" 1 !) for all n > 1. By the identity (JS]), we therefore 
conclude that the truncations and M^pp 1 ^ have the same determinant, and the 

z — 1 version of Theorem 11.11 follows. 

5.2 Refinement: proof of the MRR conjecture 

The observable t(B) for ASMs B may be included by slightly modifying the homogeneous 
limit of the IK determinant. We simply have to consider vertex weights with homogeneous 
limits in {a(r, r _1 ), b(r, r" 1 ), c(r, r -1 )} at points % = 1,2, n and j = 1, 2, n — 1 

of the square grid, and different weights {a(s, s" 1 ), b(s, s _1 ), c(s, s -1 )} for the last column 
% = 1, 2, n and j = n. Defining further 

a(s, s _1 )6(r, r _1 ) 
b(s, s _1 )a(r, r -1 ) 

gives an extra contribution z t( - B ^ to the ASM enumeration. 
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Adapting the method of enumeration described above, one finds that we simply have 
to change the definition of the last column of M^sm^ to include the z dependence. This 
in turn is obtained by modifying all columns of index j > n — 1 in the infinite matrix 
Masm (we refer the reader to [I] for the technical details). The result is the following: 

Theorem 5.1. The quantity (1 + v{z — l))ZASM(x,y, z) is the determinant of the trun- 
cation to the n first rows and columns of the modified infinite matrix M' ASM , with double 
generating function 

J m asm( u > v ) = i + i 1 \ — 

1 — uv 1 — xu — v — (y — x)uv 

{ v( z - 1 ) L + yu \ " v n-x L + v_y{y - 1) + v{xu - l) 



(y(z — 1) + x)u \ 1 — xu) \ x v + (y — vx)u 

The prefactor (1 + v{z — 1)) is ad-hoc and comes from a modification of the columns 
n and higher in the infinite matrix to make it simpler. Note that the new infinite matrix 
M' ASM has an explicit dependence on n. 

Likewise, keeping track of the observable M(A) in a DPP A is easy. The lattice path 
formulation still holds and yields a LGV-like determinant as well, but for a modified n x n 
matrix D^j, identical to Dij of (14. 3[) for i — 1,2, n and j = 1, 2, n — 1 and with a 
different last column, explicitly depending on z. 

The latter dependence is the result of decomposing further the piece (iii) of the DPP 
(see Sect J4.2p . when j = n, into {iii — a) and [Hi — b) and attaching weights as follows: 
(iii-a) the piece of the path from (£, £+ 1) to its first vertex on the y = n line at (m, n) with 
a total of ("^r^ 1 ) paths all with weight y l ~ m and (iii-b) the straight path from (m, n) to 
(0, n) with an extra weight of (yz) m+1 . This gives: 

/ k\ u , fn — m — 1 



This leads to the following: 

Theorem 5.2. The quantity (1 + z/(z — \))Zdpp{x, y, z) is the determinant of the trun- 
cation to the n first rows and columns of the modified infinite matrix M' ASM , with double 
generating function: 

JM' DPP {U,V) = + - r 

DFF 1 — zw 1 — zl — xz — w — {y — X)ZW 

_ l-v yu + u(l-xu) / 1 + " ^ 

1 — it 1 — [y{z — 1) + x)u \ 1 — xu J 

Like M' ASM , the infinite matrix M' DPP has an explicit dependence on n. 
The proof of the complete Mills- Robbins-Rumsey conjecture follows from the following 
elementary lemma, easily proved by direct computation: 

Lemma 5.3. We have the relation: 

(1 + ( y _ xv _ !)„)(! _ v)f M , ASM (u, «) = (!- u)(l + (i/ - l)v)f M , npp (u, v) 
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As explained above, such a relation between the two infinite matrices M' ASM and M' DPP 
guarantees that the determinant of their truncation to their first m rows and columns 
coincide, for any m > 1, so it holds in particular for m = n and Theorem 11.11 follows . 

5.3 More refinements 

In [5] a further observable was considered for ASMs and DPPs. For any ASM B let b(B) 
be the number of entries to the right of the unique 1 in the bottom row of B. For any 
DPP A of order n, let P(A) be the number of parts equal to n — 1 plus the number of 
rows of length n — 1. Defining the two following partition functions: 

Z^ M (x,y,z,w) = x N ^y Inv(B) - N{B) z t(B) w b ^ 

nxnASM B 

zP PP (x,y,z,w) = x s ^y NS(A) z M{A) w PiA) 

DPP A of order n 

for respectively ASMs of size n and DPPs of order n, we have: 
Theorem 5.4. |5]/ We have the identity: 

Z asm( x > V> z , w ) = Z dpp( x > V> z y w ) 

This was proved in [5] by showing that both functions Z^g M (x, y, z, w) and Z^ pp (x, y, z, w) 
obey the following relation as functions of z, w, n: 

(z-w)Z( n) (z,w)Z ( - n ~ 1 \l,l) = (z-l)wZ^\z,l)Z^ n - 1 \l,w)-(w-l)zZ^ n - 1 \z,l)Z^(l,w) 
in both cases as a consequence of the Desnanot-Jacobi identity (11. ip . 

6 Conclusion 

6.1 ASM,DPP,TSSCPP, FPL,DPL, etc. 

In these notes, we have detailed the refined enumeration of ASMs and DPPs and estab- 
lished an identity between them. 

One could think of further refinements, leading eventually to a bijection between these 
objects. This is however only the tip of a much larger iceberg (see FigJS]), which on the 
pure combinatorics side involves other objects: the so-called TSSCPPs (Totally Symmetric 
Self-Complementary Plane Partitions) which are yet another kind of plane partitions, with 
a formulation in terms of different configurations of non-intersecting lattice paths (see for 
instance [6] for a detailed account). There is also a statistical physics side, involving the 
so-called Fully Packed Loop (FPL) model on a square grid, whose configurations are in 
bijection with those of the 6V-DWBC model. The latter plays a central role in the so- 
called Razumov-Stroganov (RS) conjecture [23] proved by Cantini-Sportiello [7], relating 
its refined enumeration according to link patterns of connections of the loops around the 
grid to the asymptotic probabilities of connections of the Densely Packed Loop (DPL) 
model on a semi-infinite cylinder of finite perimeter. The latter is yet another integrable 
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Figure 3: From left to right: ASM, 6V-DWBC and FPL, all in bijection; dense loop gas (DPL): its ground- 
state/limiting probability vector satisfies the qKZ equation, the components measure FPL correlations (RS con- 
jecture), their sum matches the 6V-DWBC partition function with inhomogeneous spectral parameters Zi,u>j and 
q 3 = 1; DPPs: their refined evaluation matches that of ASMs (MRR conjecture); TSSCPPs: their refined enu- 
meration matches a sum rule for qKZ solutions at generic q and Zj = 1; Variety M 2 = 0: its degree/multidegree 
matches solutions of qKZ for q — 1. 



lattice model based on some pictorial representation of the Temperley-Lieb algebra, whose 
groundstate vector is a solution to the quantum Knizhnik-Zamolodchikov (qKZ) equation 
[9j [in] • Finally, there is an algebraic geometry side of the iceberg. For instance, the degree 
of the variety of upper triangular complex matrices with vanishing square corresponds to 
a refined enumeration of TSSCPPs, and the (equivariant cohomology) multidegree is 
obtained via a specialization of the solution to the qKZ equation jllj . 

Many of the known enumerations of the above objects involve determinants,. In a 
number of cases, these can be obtained through some application of the Lindstrom-Gessel- 
Viennot theorem. This applies to all the "free fermion" cases that are in bijection with 
non-intersecting lattice path configurations. However, both the 6V model and the DPL 
model are models of interacting fermions, in which even if there is some kind of lattice 
path formulation, the latter are no longer just non- intersecting. For instance, in the case 
of the 6V model, one can define paths going from the left border of the grid to the top 
border, by going right and up along the oriented edges as much as possible (i.e. when there 
is a choice, always go up). Such paths are now "osculating" in that two paths can bounce 
against each other at a vertex (the first going right, then up; the second going up then 
right; this corresponds to the vertex a% of the 6V model), and this configuration receives a 
different weight, interpreted as the exponential of some interaction energy. Yet, somehow, 
our formula for the refined enumeration has magically disentangled this interaction, to 
make it look like a free fermion model, via our determinant evaluation. This mechanism 
deserves to be better understood. 
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6.2 Integrabilities 



The 6V model is the archetype of 2D integrable lattice model, related to the ID quantum 
XXZ spin chain for sfa- It is know to have an infinite family of transfer matrices T(a, b, c) 
provided the Boltzmann weights a, b, c satisfy the following relation: 



a 2 _|_ 5 2 _ c 2 

A(a, b, c) = = const. 

2ab 

This constant is the anisotropy of the associated quantum spin chain. Alternatively, in 
terms of the x, y variables of (13.21) . this turns into the following "6V" variety: 

ip(x,y) = ~*~ = const. (6.1) 

as x, y > 0. 

On the other hand, the infinite matrix Masmi whose finitely truncated determinant 
gives the DWBC homogeneous 6V partition function on a grid of same size, also involves a 
transfer matrix 9 of a form analogous to that of 1+1D Lorentzian triangulations, generated 
by: 

fe(u,v) = - r — 

1 — xu — v — (y — x)uv 

The transfer matrices 9 commute for different values of the parameters x, y provided they 
belong to the following "Lorentzian" variety: 

(fix, y) = - X - = const. (6.2) 

>x 



obtained by rephrasing (12 .4p above. 

Comparing (16. ip and (16.21) . we see that (p(x,y) = ip(y,x), hence the two varieties are 
distinct! However, they do intersect. Solving for <p(x, y) = q + q~ x and ip(x, y) = p + p^ 1 
with say q,p > 1, we find that 



p(q 2 - 1) q(p 



2 



x = — r Vv 



p 2 q 2 — 1 p 2 q 2 — 1 

Conversely, any such point for p, q > 1 lies at the intersection of two "integrable" varieties 
of the form (16.11) and (16.21) . This intriguing fact deserves a better understanding. In 
particular, the 6V variety involves commutation of finite size transfer matrices, whereas 
the Lorentzian one concern matrices of infinite size. 



A Infinite matrices and truncated determinants 

Throughout these notes, we make extensive use of generating functions for infinite matri- 
ces. Let us summarize here the main definitions and properties we use. 
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A.l Infinite matrices 



We consider infinite matrices A = (cijj)ij<=z + - The very concept of an infinite matrix is 
a bit delicate to work with, for instance the product of two such matrices might not be 
well defined. This may be repaired by introducing a formal expansion parameter e, and 
associating to A the matrix A(e) = (t i+:i ai,j)i,j£Z + - The product of any two such matrices 
now makes sense in the sense of formal power series of e. Moreover, even the notion of 
eigenvector and eigenvalue make sense in this setting, provided one can show that the 
latter have formal power or Laurent series expansions in e. 

Such a construction will always be implicit (when not explicit) throughout these notes. 
For instance, the parameter g in T(g, a) of (I2.2H2.3I) plays the role of e. The "diagonaliza- 
tion" of T(g, a) along the integrable variety (p(g, a) = q + q' 1 (see Sect l2.3l) is an example 
of such extended notions of eigenvectors and eigenvalues. 

A. 2 Generating functions for infinite matrices 

For an infinite matrix A = {aij)ij ( zz + and a vector w = (u>i)i e z+ we define the formal 
generating functions 

f A {u,v) = ^ dijU'v 3 f w (u) = ^ W i U% 



with the following properties for matrices A, B and a vector w. 

1 



fi(u,v) = 

1 — uv 

f A *(u,v) = /a(v,«) 
fAw{u) = j> 1 ^- t f A {u,r 1 )f w {t) 

/dt 
-^-^f A (u,t~ l )f B (t,v) 



where the contour integral picks the constant term in t. 

We consider the lower and upper triangular matrices L(a, (3) and U(a, (3) with gener- 
ating functions 

+ H fu ^ )M= l-Pv(l + au) (AA) 

with U(a, (3) = L(a, f3) 1 . Let us also introduce the shift matrix S = (^j+i)ijGZ + and the 
transfer matrix T(a,^,j) generated respectively by: 

u 1 

fs(u,v) = - fT(a,^){u,v) = (A.2) 

1 — uv 1 — au — pv — 'juv 

In particular, we have 

L(a,0) = T(p,O,a0) U(a, 0) = T(0, 0,a/3) I = T(0,0,1) 
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The case of the infinite transfer matrix T(g,a) for Lorentzian triangulations corresponds 
to the identification: 

T(g, a) = T(ga, ga, g 2 (l - a 2 ) 

We have the following properties easily derived by contour integrals for the correspond- 
ing generating functions: 

L(a,P)L(a',P') = L (^^(l + «/?')) (A.3) 

U( a ,0)U{a>,F) = U (j^p^( l + a 'P)) (A-4) 

Ha ,0)-* = L(-l,-±) UMr = u(-y^ (A.5) 

L{a,(3)U{a',(3') = T ((3, (3', (3(3' (aa' — 1)) (A.6) 

TT( ' R'\T( R\ 1 T ( a '^' af3f3 ' (A ^ 

U («>V L W) = i^W' T {i^w^T^W-i^W') (A - 7) 

rp( a ^rpf. Of /\ _i rp ( « + V P' + l'P lj ~ (3' \ , , 

T(a,/3, 7 )T(a,/3 ,7) = — T I — , — , — - (A.8) 



1 - (3a' V. 1 - Pa ' 1 - /3a ' 1 - /3a' 

r(«,/j, 7 )-' = ^rf-i,^!) (A.9) 

a/3 + 7 V 7 7 7/ 
These hold whenever the denominators are non-vanishing. 

A.3 Commuting families and addition formulas 

Using the formulas above, it is easy to derive the following: 

Theorem A.l. The following family {T Stt (a)} ae c of infinite matrices commute among 
themselves for any fixed values of s and t: 

T s ,t{a) = T(a, sa, 1 — ta) 

We also have the following "addition" formula: 

r - I - l 1 \ f-r-i / /\ 1 ^ ( a + a' — taa' 
T s Aa)T s Aa') = - -T ' 



1 — saa' ' \ 1 — saa' 
For s = 0, we obtain a family of commuting lower triangular matrices 

L t (a) = T 0;t (a) = L Q- - i, a 

with the "addition" formula: 

L t (a)L t (a') = L t (a + a — taa) 

Changing variables from a to a, with a = 1 ~ 6 f - , and writing £ t (a) = L t (a) we finally get 
the addition formula: 

£ t (a)e t (a')=e t (a + a') (A. 10) 
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from which we deduce that £t(a) is an infinite matrix exponential. More precisely, let M t 
be the infinite matrix generated by 



(tv - l)u 
(1 — uv) 2 



Mt 



-1 



V : 





t 

-2 2t 

-3 



3t 
••••/ 



then we have £t(a) = exp(— aM t ), which by triangularity holds for any finite truncation 
as well. A similar analysis holds for T 8> t(ct), but only for the infinite matrix. Assuming 
that t 2 — As > 0, and introducing another parameter r = -i/l — As/t 2 , the relevant change 
of variables is: 

2(e rta - 1) n TM 

T ri t(a) = T Sit {a) 



a 



t(te rta (r + 1) + r - 1) 



in terms of which 



T r ,t( a ) T r,t( a ') 



1 



1 - 



(l-r 2 )(e rto -l)(e rto '-l) 
(te rta (r+l)+r-l)(te rta ' (r+l)+r-l) 



T r Aa + a 



This reduces to (1A.10I) when r = 1 (corresponding to s — 0). 



A. 4 Truncated determinants 

For any infinite matrix A = {di^)ij&z+i we denote by A^ ),n ~ 1 ^ the finite nxn truncation of 
A to its n first rows and columns, namely the matrix with entries: A' '" -1 ' = (ctij)jj'e[o,n-i]- 

In general, the matrix product does not respect truncation. However, if L, U are re- 
spectively lower and upper triangular infinite matrices, then (L[/)' 0,ri-1 l = L^°' n ~ 1 ^U^°' n ^ 1 \ 
Note that this does not hold for (UL) (see the example below). 

Let us now examine truncated determinants, namely the determinant of such finitely 
truncated matrices. By triangularity it is immediate to compute: 

det (L(a,(3)^) = (a^ k+1 ^ 2 = det (U{a, (3)™) (A.ll) 

and by the above property we deduce from (1A.6j) and (1A.11I) that: 

det [T^ k \[5^\[5(5'{aa' - 1))) = det (L(a, [5) m U{a\ = {a(3a' 73') fc(fe+1)/2 

and more generally 



det (T^ (a, p, 7)) = {a/3 + 7 ) fc ( fc+1 V 2 



whereas 



det l(U(a', f3')L(a,/3)) m 



uei ^.i yx-ppn 1-/3/3" 1-/3/3' 

(1 - (3f3') k+l 
(a/W) fc{fc+1)/2 



/U+l 
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by use of (IA.7I) . The discrepancy with the LU result is because the matrix product now 
involves all the elements of the (k + l)x infinite and infinite x(/c + l) rectangular matrices 
of the truncated product. 

The main property allowing for proving truncated determinant identities from relations 
between generating functions is the following: 

Lemma A. 2. Let L, U, A be respectively lower triangular, upper triangular and arbitrary 
infinite matrices, and let M = LAV . Then: 

M [0,k] _ jjp,k] A [0,k] jj[0,k] 

Assuming further that both L and U are unitriangular, we then deduce that det(M' 0,fe ') = 
det( J 4' 0,A: l) for all k > 0. So we will have identity between all truncated determinants of 
two infinite matrices M and A if there is a relation M = LAU for L, U lower and upper 
unitriangular infinite matrices. 
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